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ABSTRACT 

We use a compilation of cosmic microwave anisotropy data (including the recent VSA, 
CBI and Archeops results) , supplemented with an additional constraint on the expan- 
sion rate, to directly constrain the parameters of slow-roll inflation models. We find 
good agreement with other papers concerning the cosmological parameters, and dis- 
play constraints on the power spectrum amplitude from inflation and the flrst two 
slow-roll parameters, finding in particular that ei < 0.057. The technique we use for 
parametrizing inflationary spectra may become essential once the data quality im- 
proves significantly. 

Key words: cosmology: theory — cosmic microwave background 



1 INTRODUCTION 



Recent measurements of the cosmic microwave background 
(CMB) show a flat portion at low multipole number I and a 
sharp peak around I ~ 200, as well as tentative evidence for 
a peak structure beyond I = 200. This represents a tremen- 
dous success for the simplest models of the universe de- 
scribed by a flat Friedmann-Robertson-Walker metric with 
adiabatic perturbations, which are in excellent qualitative 
agreement with these observations. The power of the CMB 
is that it can be used to constrain cosmological parame- 
ters, as well as allowing us to test our assumptions about 
the form of the initial irregularities in qualitative and now 
quantitative ways. The most popular assumption concern- 
ing the initial irregularities is that they originated during a 
period of cosmological inflation (see Liddle & Lyth 2000 for 
an extensive account of inflationary cosmology). 

There have now been several papers which have 
searched for possible effects in this data from quite com- 
plicated inflationary models. One example is the inclusion 
of extra isocurvature degrees of freedom in the primor- 
dial power spectrum in addition to a dominant adiabatic 
component (Trotta, Riazuelo & Durrer 2001; Amendola et 
al. 2002), with the conclusion that the current data set is 
consistent with a sub-dominant isocurvature component (or 
even a dominant one on large scales in the case where the 
isocurvature perturbations are correlated with the adiabatic 
ones) and that the allowed values and ranges of the cos- 
mological parameters are sensitive to the type of perturba- 
tions under consideration. Another example is attempts to 
flt inflation-motivated 'features in the power spectrum' to 



the data (Grifliths, Silk & Zaroubi 2001; Barriga et al. 2001; 
Adams, Cresswell & Easther 2001). These scalar power spec- 
tra have the intrinsic property of introducing extra degrees 
of freedom — the shape parameters associated with the fea- 
ture — which can be used alter the peak heights at will. 
However it is necessary to choose the features to coincide 
with characteristic scales in the CMB power spectrum, such 
as the flrst or second peaks [early work in this direction 
(Adams, Ross & Sarkar 1997) was also motivated by possi- 
ble features in the matter power spectrum] . The CMB spec- 
trum alone offers no evidence for any extra features and so 
smooth power spectra are currently best motivated. 

It is surprising that relatively little attention has been 
paid to the CMB spectrum resulting from slow-roll inflation, 
even though these models have been the most intensively 
studied since its conception. Slow-roll inflation has acted as 
a guiding principle for inflation model builders, and is the 
simplest assumption and capable of giving excellent agree- 
ment with observations. The reason why speciflc studies of 
slow-roll inflation have been lacking is the usual assumption 
that inflation predicts a nearly power-law shaped spectrum, 
and hence that any information about inflation can be ex- 
tracted as some linear combination of the constraints on the 
two key parameters, the scalar spectral index ns, and the 
tensor fraction R, an approach used by Kinney, Melchiorri & 
Riotto (2001) and by Hannestad et al. (2002) to discuss con- 
straints on inflation. Hansen & Kunz (2001) also included 
the running of the spectral index, translating constraints on 
these parameters to place bounds on derivatives of the in- 
flaton potential. Other parameter analyses (Wang, Tegmark 
& Zaldarriaga 2002; Percival et al. 2002) have tended to fo- 
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cus on results for other cosmological parameters such as the 
densities of the various matter components, and have been 
content to use this simple parametrization. 

However, this approach ignores the fact that current 
data are only weakly constraining, and the current data set 
permits parameter regions where a significant deviation from 
a Harrison-Zel'dovich spectrum is allowed, and where the 
use of the full slow-roll power spectra is required to obtain 
robust results. The principal aim of this paper is to make 
the first direct estimation of slow-roll inflation parameters 
from CMB data. While the full predictions are relevant only 
in extreme regions of parameter space given current data, as 
the global cosmological data set improves (and in particular 
as significant observational weight develops on the high I 
part of the CMB spectrum) it may well be that these types of 
corrections take on increasing importance, depending which 
(if any) inflation models prove capable of fltting the data. 

We do not expect any dramatic new results from this 
analysis given the quality of present data, though it is a use- 
ful test of the robustness of results under more general forms 
of the initial power spectra. However it is an important test 
of principle to bring these methods to bear on cosmological 
parameter estimation, as the high-quality data of coming 
years may well require the high-accuracy description of the 
power spectrum that this approach allows. 



2 INFLATIONARY PARAMETERS 

In order to implement the inflationary cosmology into a pa- 
rameter search method, we need an adequate parametriza- 
tion of the scalar and tensor power spectra that give rise to 
the observed anisotropics. We use the recently introduced 
horizon-flow parameters (Schwarz, Terrero-Escalante & Gar- 
cia 2001) which are based around the Hubble parameter 
during inflation and its derivatives. These parameters en- 
ter directly into both the Friedmann equation and into the 
mode equations for scalar and tensor perturbations. They 
often simplify the results from analytical calculations of the 
power spectrum because of their simple definition 

eo = Hi^f/H- (1) 



where il/inf is the Hubble parameter at some chosen time 
and N is the number of e-foldings of inflation. More impor- 
tantly from the observational side, this formalism provides 
a consistent and detailed description of the shapes of the 
scalar and tensor power spectra as well as their absolute 
and relative normalizations, independent of other cosmolog- 
ical parameters (such as the cosmological constant physical 
density uja)- 

Our inflationary parameter set consists of the parame- 
ters {Hinf, ei, £2, £3, . . .}, evaluated at some particular time 
during inflation. As the amplitude of inflationary pertur- 
bations is primarily determined around horizon crossing, we 
can relate time to scale by asking when a given scale equalled 
the Hubble radius during inflation, k — aH , and we can then 
think of these parameters as a function of scale. The scale 
at which they will be evaluated is in a sense arbitrary but is 
of course most wisely chosen to be around the centre of the 
scales actually probed by the observations. This parameter 



set replaces that made up of astrophysical quantities {Viz, 
Rio{flA, h), ns — 1, dns/dln k, ut, driT /dink, . . .} that are 
used if the power spectra are taken as the starting point. 
Here Viz is the amplitude of scalar perturbations specified 
at the scale k = 0.05 Mpc~^, _Rio is the ratio of the tensor 
and scalar d curves evaluated at I — 10, and ns — 1 and 
riT are the slopes of the initial scalar and tensor power-law 
spectra. Use of the infiationary parameters automatically en- 
forces the infiationary consistency relations between scalar 
and tensor perturbations, and so there are fewer parameters 
to be fit; a full treatment of observations would also test 
these conditions though present data is not good enough for 
a meaningful confirmation. 

The current data set has very limited abilities to con- 
strain tensor perturbations and usually only R is used to 
describe them (perhaps coupled with a consistency rela- 
tion to fix the tensor spectral index). It is also common, 
though not universal, for parameter searches to neglect 
scale-dependence of the scalar spectral index. 

The power spectrum V from slow-roll infiation can be 
obtained as an expansion of the power spectrum (or some 
function of the power spectrum) in terms of the logarith- 
mic wavenumber. While the usual power-law expression is 
mostly simply related to an expansion oilnV, since it is the 
power spectrum itself which is constrained by observations 
it is most direct to expand 'P(fc) itself (Martin, Riazuelo & 
Schwarz 2000; Leach et al. 2002). The spectra of curvature 
perturbations Vtz and of tensor perturbations Vh are writ- 
ten as 

p7^,hOcfeo + blln(A) +^ln2 (^A) +..., (3) 

The detailed predictions for the bi have been calculated to 
second order in the slow-roll parameters (Stewart & Gong 
2001; Leach et al. 2002) to be 

+ - C + ^ - 7) £162 + (iC' + ^~ljel 



+ {~^C^ + ^)^2e, (4) 

&S1 = -2ei-e2 + 2(2C4-l)e? + (2C-l)eie2 

+Cel - C€2€3 (5) 
bs2 = 4ei 4- 2eie2 + £2 — £2 63 (6) 

for the scalars, and 

6to = 1 - 2 (C + 1) ei + (20^ + 2C + 4 - 5) e? 

+ (-C'-2C+f^-2)ei62 (7) 
&T1 = -2ei + 2(2C + l)e?-2(C-(-l)ei£2 (8) 

&T2 = 4e? - 26162 (9) 

for the tensors, where C = 7E+ln2 — 2~ —0.7296 is a 
numerical constant. The full parametrization of infiation as 



given here, where we truncate at 63, should remain suffi- 
ciently accurate for some time to come, quite likely including 
Planck satellite results. 

Concerning the relative normalization of tensor and 
scalars, the constant of proportionality for Eq. (^ is given 
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by the slow-roll amplitudes 



(10) 



Vnoik*) = (11) 

7reim|,j 

In the slow-roll limit the tensor to scalar ratio is given by 
Vh/V-R. — 16ei, and the more general expression can be read 
from the above expansion coefficients as 



^ = 16ei[l-f Cea 



(12) 



This logarithmic expansion of the power spectra is ac- 
curate over a wide range of the inflationary parameter space. 
However, unlike the expansion of InP(fc) (the power-law ex- 
pansion) , the logarithmic expansion of 7'(fc) can become neg- 
ative at large |ln(fc/fc«)| for large e^. This pathological be- 
haviour serves as a warning against using either expansion 
alone for robustly extracting any inflationary signal in this 
regime without cross-checks. In any case, this is behaviour 
does not occur in our analysis due to the limited range of 
scales currently probed by the CMB. The methods that we 
describe in this paper therefore complement the traditional 
approach, and a careful analysis should utilize both (Leach 
et al. 2002). 

To lowest order in slow-roll these predictions reduce to 
the well-known linearized form 



ns 



R: 



Vtz 



-261-62, 

-2£i, 



= 16ei. 



(13) 
(14) 

(15) 



Written is this way it is clear that an observational strat- 
egy would ideally use information from the tensor sector 
(primarily the tensor amplitude) to break the degeneracy 
between ei and £2 for the scalar spectral index. The ef- 
fect of increasing ei is to boost the large-angle anisotropics 
via the tensor amplitude while simultaneously tilting down- 
wards the small-angle anisotropies via the scalar tilt, which 
act roughly constructively. 

In slow-roll inflation it is possible to relate the horizon- 
flow parameters to the shape of the inflationary potential 

~ J^V, (16) 



ei 



£2 



3m| 

■mil 
16tt 
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4-K 



327r2 













V"' 
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V"'V' 


-3^ 


1/2 


V 



(17) 



(18) 



(19) 



3 OBSERVATIONAL CONSTRAINTS 

We fit to a CMB dataset comprising data from COBE, 
BOOMERanG, Maxima, DASI, VSA, CBI and Archeops. 
We follow the method of Lesgourgues & Liddle (2001) in 
defining a simple with penalty functions over De = 
{ATif = ^(^-f l)C£/27r for the COBE (Bennett et al. 1996; 



Tegmark & Hamilton 1997), BOOMERanG (Netterfield et 
al. 2002), Maxima (Lee et al. 2001) and Archeops (Benoit 
et al. 2002) data of the form 



E 



- (1 + car + &o-b,i)Df 



(20) 



where for COBE we have o-b,i = (Jc = 0, BOOMERanG = 
0.20 and o-b,i = 0.43 x 10"'^'^^ Maxima crb,i = 2 x IQ-^i^'^ 
and (Jc = 0.08, and Archeops CTc = 0.14 and CTb.i ~ 0. For 
the DASI (Halverson et al. 2002), VSA (Scott et al. 2002) 
and CBI (Pearson et al. 2002) ('mosaic' configuration, odd 
binning up to ^ = 1500) data we define the to be 
2 



X^ = {Dl^°° - {l + c(Jr)D°^'') M.T^ 



(l + cac)I>rH-c-^ 



(21) 



where ac = 0.08,0.06,0.10 respectively and the correlation 
matrix is given by 



Mij = AA'KjADj + alDiDy 



(22) 



We assume gaussian window functions, except for the CBI 
data for which the published window functions are used, and 
the Archeops data for which the window functions are taken 
to be top-hat functions. All error bars are approximated 
as symmetric, taking the larger error bar when they are 
not. We analytically determine the coefficients bboom, Cboom, 
&max, Cinax, Cdasi (simply by simultaueously solving equations 
such as dx^ / dbhoom = 0) and then sum the x? over all the 
experiments. 

The CMB models were calculated using the January 
2002 version of CAMB (Lewis, Challinor & Lasenby 2000) 
which can be easily be modified to incorporate absolute and 
relative normalizations of the scalar and tensor power spec- 
traj^ We note in passing that indeed the correct amplitude 
in the tensors is more important that the correct tensor tilt. 
Even when we move the pivot scale away from COBE scales 
where the tensor spectrum is physically relevant, the tensor 
to scalar ratio only runs weakly back to the COBE scales 

^ = 16ei[l + (ln(fc/fc.) + C)e2 + ...]. (23) 

The actual shape of the tensor spectrum will have very little 
effect at this stage, and so we are certainly justified in taking 
the same pivot scale for both the scalar and tensor spectra. 

We examine a set of spatially-flat cosmologies using the 
parameters {ojb, wm, ti^A, e~^'^, Ptj, ei, £2} where Ui = Qih^ 
measures the physical matter density and r is the optical 
depth to the last-scattering surface. Our grid runs overs the 
range 0.006 < t^B < 0.045, 0.03 < t<;M < 0.28, 0.0 < uja < 
1.2, 0.28 < e-'^ < 1.0, 7 < Vtz X lO"' < 48, 0.0001 < 
ei < 0.07, —0.4 < £2 < 0.3, with a uniform spacing of 
8^ X 7 X 22 X 8 X 11. It is worth emphasising that we have 
included Vtz as a full parameter instead of using an analytic 
procedure to choose the best amplitude from the data and 
each model. Apart from our obvious interest in primordial 
physics, Vtzc"^'^ is a parameter that will be determined with 
high accuracy in the future and so should be included in the 
analysis. 

^ A module to directly input the predictions of slow-roll in- 
flation to the CAMB program is available to download at 
www. astronomy . Sussex . ac.uk/~sleach/inflation/ 
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The dependence of the power spectrum on the third 
slow-roU parameter £3 is very weak, given the current data 
set, and so we fix ea = while stiU using the second-order 
expressions, which give a better representation of the power 
spectrum at large ei,e2 and |ln(A:/A:«)| than the first-order 
expressions. Our decision to truncate the slow-roll expan- 
sion after e2 is based on subtle considerations. As the series 
expansion is an infinite one, it clearly has to be truncated 
somewhere; if we did include ea the same question of whether 
to include or not would then arise for £4. We believe that a 
suitable criterion is that one ought not to include parame- 
ters which the data cannot distinguish from zero, i.e. we im- 
pose a null hypothesis that those parameters are zero. While 
that criterion is driven by data alone, it makes sense in the 
context of inflation because one does expect the parameters 
to be small, and if the data is unable to constrain them 
to be small, then the parameter space explored would in- 
clude models which would not be satisfactory inflation mod- 
els. Unfortunately this criterion cannot be rigidly applied, 
as even ei and £2 cannot presently be distinguished from 
zero; thirty years of large-scale structure studies have not 
excluded the Harrison-Zel'dovich spectrum. However if in- 
flation is to be meaningfully tested, eventually at least one of 
these parameters needs to be detected with a non-zero value 
and so we do include them in our current treatment. This is 
appropriate as inflation models do predict values of ei and 
£2 comparable to constraints from present data. By contrast, 
typical inflation models predict £3 to be much smaller than 
the present limits from data (see e.g. Hansen & Kunz 2001). 
Whether experiments such as Planck can detect a non-zero 
£3 remains to be seen. 

We fix the pivot point to fc* = 0.05Mpc~^, ignoring the 

1/2 

cosmological dependence, proportional to lu^ , of the range 
of scales probed by the CMB. The Hubble parameter is an 
auxiliary parameter given by 



h = V UJM + OJA 



(24) 



Throughout this paper we examine the effect of applying the 
HST Hubble parameter prior of ft = 0.72 ±0.08 at la, which 
acts as a constraint on the total physical matter density for 
our fiat cosmologies. 

When presenting parameter constraints, we simply min- 
imize the over unwanted parameters. In principle one 
should consider performing a cubic spline over the of 
our coarse grid of models in order to examine the shape of 
the likelihood function e ^ ' and to subsequently perform 
a marginalization procedure (Efstathiou et al. 1999). How- 
ever, minimizing the reproduces the basic shape of the 
allowed region in parameter space (Tegmark & Zaldarriaga 
2000). Our best-fit model has = 51.1 and so we plot the 
2(7 and Scr contours for a x^ distribution with 51.1 degrees 
of freedom. We omit the la contour, as little useful infor- 
mation is conveyed by this level of certainty in the context 
of CMB parameter searches. We regard the models enclosed 
by the 2a contours as representing good fits to the data, and 
the 3cr contours mark out the region where tension between 
the data and the models begins. 

Turning firstly to three of the standard cosmological pa- 
rameters, the densities of the three main components, we see 
in Fig. |l| impressive agreement with the standard BBN value 
of ujs = 0.020 ± 0.002 at 95 per cent confidence (Buries, 
NoUett & Turner 2001). This acts as a useful consistency 
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Figure 1. Three views of the constraints on wb, i^m and oja. The 
dashed lines show the constraints from the CMB alone, while the 
solid lines include the HST prior on the Hubble constant. 



check on the assumption of adiabaticity, given that the in- 
clusion of a subdominant isocurvature mode tends to widen 
the allowed range in ujs, (Trotta et al. 2001). For the cur- 
rent dataset we observe a weak correlation between the con- 
straints on ul!m and uja- The main effect of applying a prior 
on the Hubble parameter, h, is to rule out models with a 
large physical density in cja, and to close the contours near 
UJA = 0. When quoting parameter constraints we include the 
HST prior. 
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Figure 2. Constraints on the amplitude V-n and the first slow- 
roll parameter ei, the solid lines again including the HST prior. 



We now turn to the inflationary parameters, which are 
the main focus of our study. In Fig. ^ we plot the constraints 
on Vti and ei. The general trend is that as we increase the 
tensor component, R — 16ei, which contributes at low £, 
then we must decrease the scalar normalization in order 
to continue to fit the low £ data. Including Vn as a full 
parameter in the analysis acts as a useful warning against 
including higher-order parameters in the power spectrum; 
until we have pinned down the amplitude of perturbations 
at some given scale, there is little sense in trying to mea- 
sure the higher derivatives of the power spectrum such as 
the running of the spectral index etc, under the assumption 
that these higher derivatives are weak. Reading off the 2-a 
bounds on Viz and ei we find 



IKVtzX 10'" < 42, 
ei < 0.057. 



(25) 
(26) 



These are the main results of this paper, and are in good 
agreement with the analysis of Lewis & Bridle (2002). The 
upper bound on ei is consistent with the inflationary hy- 
pothesis, which requires ei < 1 for inflation to occur. It isn't 
possible push this limit much further down using CMB data 
from ground and balloon observations due to the calibration 
uncertainties in the high I region which mask the primary 
effect of 61 , which is to boost the low Ces relative to the high 
Cis. 

Reading off the values of Viz and ei at the tip of the 
2-0" contours of Fig. ^ allows us to place a constraint on the 
energy scale of inflation 

— < 1.5 X 10-^ (27) 



mpi 



V 



1/4 



< 2.9 X 10"'' GeV. 



(28) 



Knox & Song (2002) and Kesden, Cooray & Kamionkowski 
(2002) have calculated a limit on V^l/f* below which tensors 
can not be detected directly in the B-mode polarization of 
the CMB due to a contaminating signal from lensing of the 
i?-mode along the line of sight. Combining this with our 
upper limit then we find that the inflationary energy scale 
must lie in the range 



3 < 



V 



1/4 
nf 




Figure 3. Constraints on ei and 62- The data are consistent with 
a scale-invariant, tensorless spectrum (ei << 1, |e2| << 1), but a 
significant tensor fraction and tilt are still permitted by the data. 
The dashed line, ei = —62/2, indicates inflationary models with a 
scale-invariant power spectrum. To the right the spectra are red, 
to the left they are blue. 




Figure 4. Constraints on Vtz and the optical depth r. The thick 
dashed curve is given by IO^^'Ptjc"-^'^ = 15. 



in order to directly detect the tensor spectrum via the B- 
mode polarization of the CMB. 

In Fig. ^ we plot the constraints on ei and £2, the flrst 
two slow-roll parameters, showing that a wide range of slow- 
roll inflation models fit the present data. The data are con- 
sistent with a scale-invariant scalar spectrum with no tensors 
(ei <^ 1, |e2| ^ 1). We can also read off an approximate and 
loose bound on the second slow-roll parameter 



0.31 < e2 < 0.2. 



(30) 



IQiSGeV 



< 29 



(29) 



To guide the eye we plot the line ei = —£2/2 along which the 
inflationary scalar power spectrum is approximately scale- 
invariant. As can be expected, the contours lean in the same 
direction as this line reflecting the main inflationary degen- 
eracy of Eq. (|l|). 

Finally, in Fig. ^ we plot the constraints on Vn and op- 
tical depth T revealing the anticipated degeneracy between 
these two parameters, both of which affect the amplitude of 
the acoustic peaks. We derive the limit 0.4 < e~^^(< 1.0) 
corresponding to an upper limit r < 0.45. 
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4 SUMMARY 

We have implemented the detailed second-order predictions 
for the inflationary power spectra, given by equations 
to into a CMB parameter search method, using the 

logarithmic expansion of P{k), Eq. (^, for the first time. 
Although the present data set cannot hope to actually mea- 
sure thp ^»rpaV rirnninp; nf tViP appptral inHpY inHnrpH hy a 



significant tensor component, we derived sensible limits on 
the power spectrum amplitude Vnik = 0.05 Mpc~^) and the 
first two parameters, ei and 62, by assuming the running to 
be weak, which was achieved by considering models with 



ei < 0J07, -0.4 < £2 < 0.3 (and fixmg 63 = 0). We also 
derived a sensible limit tor ujh which acts as a usetul con- 
sistenc y check on the assumption of adiabatic perturbations 



with ar approximately power-law form. 



While the results of the present paper do not add much 
to existing studies parametrizing the spectra as power-laws, 
our paper represents an important point of principle in im- 
plementing precise slow-roll inflation predictions for the first 
time. As the global dataset improves, including MAP and 
then Planck data, it is quite likely that these techniques are 
required to ensure robust estimation even of cosmological 
parameters such as the densities of the various components. 
Further, these techniques will be essential to squeeze the 
maxim um possible amount of information out of the data re- 



garding infiation, should slow-roll infiation continue to give 
the simplest viable interpretation of observational data. As 
the dat a sot improvoE, it will bo intoroeting to open up tho oj. 



direct icin as well as PYplnring the possibility nf a npgligihle 
tensor prior, in order to differentiate between infiationary 
models as effectively as possible. An interesting goal would 
be to determine the signs of both e2 and ea, which would 
allow us to immediately rule out three quarters of all single- 
field slow-roll infiation models. 
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